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ABSTRACT 

Two analyses are presented to study the effects of 
restrained shrinkage on reinforced concrete members. The 
first 1S a uniaxial shrinkage analysis of a symmetrically 
reinforced, completely restrained slab element. A procedure 
for calculating induced tensile stresses and spacing of 
cracks due to restraint of shrinkage is presented. The 
Peduct1onmotax tal stit fnessiduetarcracking is) basedscnman 
assumption of linear steel stress distribution within the 
zone of influence of cracking. A parameter acu was 
conducted to investigate various factors affecting the 
member's behaviour. 

The second analysis involved the application of the 
finite element method in analyzing two-way continuous slab 
systems. Shrinkage cracks are assumed to exist in the slab 
prior to analyzing the slab for transverse loads. Methods 
for evaluating the reduced flexural stiffness due to 
shrinkage cracking and cracking due to transverse load are 
given. A brief investigation was performed to study the 
effects of number of shrinkage cracks and variation in 
effective tensile strength on slab deflections. It is 
proposed that shrinkage restraint be incorporated in a plate 
bending analysis by using a reduced effective modulus of 
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Coefficient used to determine cracked length 


inmecouem NO-Slip sapproacn. 


Coefficient used to determine cracked length 


in the 'slip' approach. 

Reduction factor to account for compressive 
steel and movement of neutral axis in 
computing long-time deflections. 

Average axial stiffness of a slab element 
after formations of one crack, two cracks 
and m cracks, respectively. 


Column stiffness. 


Axial stiffness in the cracked region of a 


Slab element. 


Axial stiffness of uncracked section of a 


slab element. 


Stiffness of a slab system with partial 


fx by 
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L Clear span length of a slab element. 
L Heuchia of. column. 


Ler Cracked length within which the stiffness 
varies from the stiffness provided by 
reinforcement alone at the crack, to the 
full stiffness of the steel-concrete 
composite member at a distance Ley/2 from 


the crack. 


Ly Development length. 

Lelip Slipping length. 

LL Live load. 

m Number of cracks formed in a slab element. 
M Moment. 

Me Cracking moment. 


Moment in x-direction (reinforcement 


parallel to x-axis) shown in Figure 3.2. 
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Twisting moment shown in Figure 3.2. 


Moment in y-direction (reinforcement 


parallel to y-axis) shown in Figure 3.2. 
Modular ratio. 

Tensile force applied to a slab element. 
Compressive force in reinforcement. 
Tensile force carried by concrete. 


Net induced tensile force ina slab element 


due to restrained shrinkage. 
Tensile force carried by steel. 


Uniformly distributed load applied to slab 


shown ine Fagure eee 
Crack spacing. 


Initial time interval used in the bisection 


method. 
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(re eter, ol Time interval obtained by the bisection 


method after k iterations. 


toe te hel Times intervalwobtained by the bisection 


method after k+1 iterations. 


T Maximum tensile load that can be carried by 


concrete. 

TTOL Acceptable tolerance in time. 

Yt Distance from centroidal axis of gross 
section, neglecting the reinforcement, to 
extreme top fibre of concrete member. 


a Axial stiffness reduction coefficient. 


ap Flexural stiffness reduction coefficient for 


precracked element. 
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ay Flexural stiffness reduction coefficient in 


x-direction. 


a y Flexural stiffness reduction coefficient in 


y-direction. 


Yyy Shearing strain. 

A Displacement due to shrinkage of a released 
support slab element. 

Jaye a= fae Displacement due to shrinkage of a released 
Support slab element with one crack and two 
cracks respectively. 

Dep Creep deflection for all sustained loads. 

Acr Elongation of cracked region under load P. 

4, Elongation of uncracked region under load P. 

Aj Immediate deflection. 
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(A;+)p Immediate deflection due to all sustained 


load. 


Acslab Elongation of a released support slab 


element with partial fixity if restraint is 


reapplied. 
At Increment of time-step. 
Bin Gp Ge Average unit strain in slab element after 


the =formatwvons or the fbirstecrack,, the 
second crack and the mth Gracky, 


respectively. 


fave Average unit strain in reinforcement. 

Ec Unite *strato®in*concrete. 

Es Unit strain in steel. 

Ech Unrestrained shrinkage strain. 

€shis €Sn2, Unrestrained shrinkage shrain at which the 
€shm first crack, the second crack and the m‘? 


crack form in a slab element, respectively. 
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Ultimate shrinkage strain. 


Normalwstrainsparallelato x-axis’. 


Normal strain parallel to y-axis. 


Pousson! Se ration 


Poisson's ratio in x-direction. 


Poisson's ratio in y-direction. 


Tension reinforcement ratio. 


Compressive reinforcement ratio. 


Steel stress at a crack. 


Steel stress at uncracked section. 


Normal stress parallel to x-axis. 


Normal stress parallel to y-axis. 


Average bond stress. 
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Shearing stress. 


Ultimate bond stress. 


Reinforcing bar diameter. 


Average curvature within precracked region. 


Curvature 


Curvature 


Shrinkage 


Curvature 


at cracked section. 


at uncracked section. 


curvature. 


in x-direction. 


Shearing curvature. 


Curvature 
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1. INTRODUCTION 


1.1 Introductory Remarks and Literature Review 


The current practice of ultimate strength design of 
reinforced concrete structures has resulted in more flexible 
and slender structures than in the past. The assessment of 
the performance of the structure at the working load levels 
thus becomes an extremely important consideration. Although 
strength requirements may be satisfied, cracking and 
deflections at the working loads may be excessive. Cracking 
may be excessive if the steel stresses are high or if the 
reinforcing bars are not properly distributed. Deflections 
may be critical when shallow sections, which are the case in 
ultimate strength design, are used and high stresses are 
present. 

The concern for serviceability requirements is 
reflected in the current ACI Building Code, where methods 
are suggested to compute both short-term and long-term 
deflections. The Code however does not explicitly consider 
erackindgsofmstructuressthateofteng occursgat anpearly age of 
construction either due to construction loading or due to 
restraint of shrinkage. Failure to recognize these effects 
may lead to unconservative predictions of two-way slab 


deflections as revealed by some investigations (Heiman, 1974 
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and Rangan, 1976). 

Observed cracking uSually occurs for the following main 
reasons (Scanlon and Murray, 1982 and RILEM Committee 42- 
CEA, 1981): 

1. Settlement of the newly placed concrete. 

2. Early shrinkage and thermal variations. 

3. Higher than design moments that occur adjacent to 
columns. 

4, Excessive construction loading due to shoring procedure 
that results in overstressing of the slab before the 


concrete can reach its design strength. 


A survey carried out by Mayer and Rusch (1967) has 
indicated that excessive slab deflection was the most common 
cause of damage to reinforced concrete structure. 
Insufficient consideration of creep and shrinkage on slab 
deflections was one of the causes cited for damage. 

The ACI Building Code recognizes shrinkage effects on 
deflections in terms of shrinkage curvature due to warping 
that arises from nonuniform shrinkage. However, in 
continuous reinforced slab systems, tensile stresses are 
induced by shrinkage of the concrete when shrinkage strains 
are reStrained. Stiff support elements provide the 
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value of the induced tensile stresses may well exceed the 
tensile strength of the concrete, particularly at the early 
stages of hardening, resulting in additional cracking of the 
reinforced slab system. The overall stiffness of the slab 
system is thus reduced, causing an additional slab 
deflection. The current ACI Building Code calculation 
procedures may, therefore, underestimate slab deflections if 
the effects discussed above are not properly considered. 

Cracking in certain regions reduces the slab stiffness 
and complicates the analysis of reinforced concrete slab 
systems. Classical plate-bending solutions cannot readily 
deal with slabs with varying stiffness. The introduction of 
the finite element method solves this problem by dividing 
the slab into finite regions allowing the stiffness to be 
varied from one element to another. 

Successful modelling of cracking behaviour in slab 
systems has been made by Jofriet and McNeice (1971). A 
bilinear moment-curvature relationship formed the basis of 
the model of a cracked region. Cracks were assumed to 
iMmittatesat a direction perpendicular tomthe major principal 
moment and unalter during any increase of load. Procedures 
were developed to give equivalent steel areas normal to the 
Clack smrOrmc rect substi cublOonminet nescalculatlonmore une 


effective moment of inertia of a cracked region. 
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A much simplified method was developed by Scanlon and 
Murray (1982), who included the effects of cracking by 
reducing the values of the elastic constants. Moments in 
the x and y directions were checked. If the cracking moment 
in either direction was exceeded, reduction in stiffness was 
made uSing the ACI Building Code equation for effective 
moment of inertia. The procedure thus assumed that cracks 
were oriented parallel to and perpendicular to the direction 
of the reinforcement. The reinforcement was placed in the x 
and y directions of the global co-ordinate system so that 
transformation of steel areas was not necessary. 

The additional shrinkage cracks that may occur in slab 
systems were not modelled in the two forementioned 
procedures. A finite element model that includes both 
cracking due to transverse loadings and due to shrinkage 


restraint is, therefore, required. 


1.2 Objectives 


The objectives of this investigation are: 

1. To investigate the effects of shrinkage cracks on 
flexural stiffness of slabs, 

2. to analyze reinforced concrete slabs that are subjected 
to shrinkage cracking using the finite element method, 


and 
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3. to assess a simplified procedure for calculating 
deflections which includes.the effects of restraint 


Stresses due to shrinkage of concrete. 


1.3 Outline of Content 


In Chapter 2, the analysis procedure of a member 
subjected to uniaxial shrinkage strain is described. The 
length and uniaxial stiffness of a cracked region is 
modelled. The induced tensile stresses in the concrete due 
to restraint of shrinkage are evaluated. Checking the value 
of the induced tensile stresses against the tensile strength 
of concrete, the number of cracks developed and the 
shrinkage strains at which cracks formed are also 
determined. A parameter study is then carried out to 
investigate the different effects of variations in several 
Parameters. 

The finite element procedure of the analysis of a slab 
system subjected to bending is described in Chapter 3. 
Cracking due to transverse loading and due to restraint of 
shrinkage is modelled. To complete the analysis, additional 
deflections due to shrinkage curvature and creep are 
approximated. The results of the analyses of three 
reinforced concrete slabs are used to assess a simplified 


method for evaluating slab deflections subjected to the 
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effects of early shrinkage cracking. 
A brief summary and recommendations for further 


Besearchearesinc muded sim Chapter 4. 


2. UNIAXIAL SHRINKAGE ANALYSIS 


2.1 Introduction 


Stresses are induced in concrete when volume change due 
to shrinkage is restrained. Restraint may be provided by 
Support conditions and/or by reinforcing bars. Tensile 
Stresses developed may exceed the tensile strength of the 
concrete, particularly at the early stages of hardening. 
When the tensile strength is exceeded, cracking occurs. 

In this chapter a symmetrically reinforced slab element 
fixed at both ends and subjected EO.UNLCOLMM weuntaxdaL 
Shrinkage is analyzed. For a given shrinkage strain versus 
time history, the effect of progressive reduction in element 
stiffness as cracks form is taken into account in 
determining the spacing of cracks in the slab element due to 
restraint of shrinkage. Results obtained are used as a 
basis for specifying the extent of cracking resulting from 
restrained shrinkage in a Slab prior to analyzing the slab 


for transverse loads. 


2.2 Assumptions and Limitations of Analysis 


The analysis has the following assumptions and 
limitations: 


1. All concrete properties except shrinkage strain are 
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assumed to be constant with time. An example of a 
shrinkage strain versus time curve is shown in 
bel.gur ees2unlr. 

2. A crack forms when the tensile stress in concrete 
exceeds the specified tensile strength. 

3. Creep effects are not included explicitly in the 
analysis. 


4, No yielding occurs in reinforcing bars. 


2.3 Shrinkage Analysis 


2.3.1 Outline of Analysis 


In the following sections, the shrinkage analysis 
procedure of a symmetrically reinforced slab element with 
full fixity at ends will be presented. The slab element is 
first considered both at the uncracked state and at the 
first and second cracking states, then the analysis is 
generalized to consider the formation of any number of 


cracks. 


2.3.1.1 Uncracked state 


Consider the member in Figure 2.2a. Release one 


support and allow the concrete to shrink freely by an amount 


A = EshL 
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where esh is the shrinkage strain that occurs in concrete. 
In order to force the steel bar to shorten by the same 


amount, a compressive force 
Omen sEcesh 


must be applied to the bar as shown in Figure 2.2b. The 
concrete 1S now unstressed and the bar is under compression 
when one support is released. 

For compatibility, a force P must be applied to bring 
the member back to its original length ( Figure 2.2c ). The 
applied load P is carried partly by the concrete (Pc) and 


partly by the steel bar (P,) at any cross section. 


Hence, he te ue te 


= yet a tryslarsis/s 


Before cracking occurs, there is no slipping between 


concrete and steel. 


Hence, Esh 
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Therefore, P = AcE¢ Esh ou As Es Esh 
= (l+np)AcEcésh 


= AeEcésh 


where Ae = (1t+np)Ac 
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A, is the equivalent transformed area of the member if 
the steel area A, is replaced by an equivalent amount of 
concrete area npA.. 


The net force in the member is 


Pret = Pp te Po 
= (ACEcésh + AsEsesn) + (-AsEsesn) 


= Vey 


Concrete stress in the member is 
fo = Ecésp = P/Ag = Pc/Ac C2 ile) 


Before cracking the net stress in the steel is zero. 
It may be noted that the same result is obtained if the bond 
between steel and concrete is released as well as the 
restraint at one of the fixed ends. In this case Py = 0 and 
the net concrete and steel stresses are aS indicated above 


when compatibility of deformations is applied. 


2.3.1.2 Formation of the first crack 


The concrete stress continues to increase as shrinkage 
Continues suntils the concrete tensilesstrengtn t+e1s reached; 
Causingmrormation of the first) crack. | Thesshrinkage: Strain 
atewhichithesfirstecrack forms .(91.6e-whengt, =f) jatcom 
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Cracking causes a reduction in member stiffness with a 
resulting decrease in stress in the member. Figure 2.3 
shows the variation in concrete tensile stress with 
increasing shrinkage strain’. slnitaaliy>eini the: uncracked 
state, the tensile stress increases linearly with shrinkage 
Strain along line OA. At point A the tensile strength is 
reached and the first crack forms. The reduction in 
stiffness causes a decrease in stress to fc; ( point B ). 
As shrinkage continues to take place, the concrete stress 
increases again along line BC up to point C where the 
tensile strength is again reached and a second crack is 
formed. This sequence continues until the shrinkage strain 
reaches its ultimate value. 

The residual concrete stress fc; can be obtained using 
the following procedure. 

Se Deel. Figure 2.4a shows the member with one crack 
formed. As a result of the crack, the axial 
stiffness is reduced over a length Lcr. Within 
this length, the stiffness varies from the 
stiffness provided by steel alone at the crack, 
to the full stiffness of the steel-concrete 
composite member at a distance Lcy/2 from the 
Crack. 

If slipping between concrete and steel 
occurs within some distance Ls]jp due to bond 
failure around the first crack, then shrinkage 


within this slipping region will not induce 
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Significant stresses. Hence, only shrinkage in 
Phegnongsiip region (lic nets considered. 
Consider again the released support member 
with one crack formed as shown in Figure 2.4b 
and allow the concrete to shorten freely by the 


amount 
A, = esh;(L-Lslip) 
Define the average member strain to be 
Che esa mic) Gare 


Therefore A, = €,L 


BhenstheRcompressivemroOrces > nequined ato 


shorten the reinforcing bar by the same amount 


zs 
Po = = AsE<sA,/L 
SS AsEseé 
Sto pe. A force P is applied to bring the member length 


back to L. Because of the reduction in 
stiffness over Ley the displacement A, is 
considered in two parts. The cracked region 


elongates by an amount 


Acr = Plcr/Kcr 
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where Kcy iS the reduced stiffness over the 
length Lcy. 
The gross section also elongates over a 


Teng thot (aiae— lca) aDyeaneanount 
Ag = P(L-Ler)/Kg 
where el EVAN 5 
Meom compatibility, 
Bo = Acr + Ag 
Therefore, PL/K, = PLer/Kcr + P(D-Low)/ Ka 
in which K, is the overall member stiffness with 
one crack. 
Rearranging the above expression, 


K = L 
(L-Ley)/Kg + Ler/Ker 


Similarly, for a member with m cracks, 


Ki tee t= SG he Pee, OE (2 4) 
(L-mLey)/Kg + MLcoy/Key 


where m= number of cracks in member 
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Step 4. 
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From Equation 2.1, concrete stress in gross 


Sect onmrs 


fo: = P/Ae 
= K,€1/Ao 
Sees ia be pe PG) 


Agi 


The net force in the member is 


Phetyesn Hehe 


K,€,; si As Es €1 


(Kes AgiBsi) ey 


(24) 


After the first crack forms, the tensile stress 


in steel is no longer zero. At the crack the 


force Pnet iS carried entirely by 


producing a net tensile stress of 
09s,cr = Pnet /As 


Outside the cracked zone Ley, the 


steel is compressive and is given 


the steel, 


C257 


stress in the 


by 
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The variation of steel and concrete stresses along the 
length of the member after the first crack forms is shown in 


Figure 2.4c. 


2.3.1.3 Formation of the second crack 


Assume shrinkage continues to take place after the 
first crack forms. From Figure 2.3, when the shrinkage 
Strannmeinereasesmirommes arto feichye the cConcretesstress 
increases along line BC where the member has a stiffness of 
K,. At point C, where ech = €shz, the concrete tensile 


Strength is reached, i.e. 


From Equation 2.5, substituting eshi by e€sh2, 


fe = Kiesh2(L-Lslip) 


Aab 


Rearranging, the shrinkage strain at which the second crack 


forms is 


€sh2 = Ae ftL (293) 
K, L-Lejip 
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After the second crack is formed, the member has a 
Sr ELNesseOtak wrandmas Sl Uppinglengtn or 2Lslip- Again 
release one support to allow the member to shrink by the 


amount 


A> = €2L 


where ai esh2(L-2Lslip)/L Lae) 


The force P applied to maintain compatibility is given 


by 


8) 
| 


— Keke 


= K2€2 GO 2ealy) 


where Kz can be obtained from Equation 2.3 by letting m = 2. 
The residual concrete stress at the uncracked section 


is 


P/Ae 


K2€2/Ap (2Rnii2) 


The net force in the member is 
Poet. =) sbactero 


where Py is the compressive force required to shorten the 
reinforcing by the amount A; when the restraint is released 


at one of the fixed ends and is given by 
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= = AcsEcé2 
Therefore Pret = K2éz2 - AcEceé2 
sent K e- yA Bee (2.13) 


The steel stresses at the crack and at the gross section can 


be determined from Equations 2.7 and 2.8. 
2.3.1.4 Formation of any number of crack 
In general, the shrinkage strain egpm at which the mth 


crack forms is given by the generalized form of 


EquatuonsSec.2eands 2.9, 


Kim - 1 ST ie MENA 
where Ky-1 iS obtained from Equation 2.4 by substituting m 
egy (ne 


th crack is confirmed by checking 


Formation of the m 
that Eshm < Eshue 
The residual tensile concrete stress iS given by the 


generalized form of Equations 2.5 and 2.12, 
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The net force in the member when the crack forms is 


derived from Equations 2.6 and 2.13, 


eis Pnhet = (Km - AsEs)em e2e lath) 


The steel stresses can be determined by substituting 
the appropriate values of concrete stress and Prat into 


Bouationsec. sand 260. 
2.3.2 Shrinkage Crack Model 
2.3.2.1 Reduced Stiffness of Cracked Region 


In the above analysis procedure, the development of 
progressive cracking as shrinkage takes place depends on the 
effective stiffness Keve over the region Myer WR Eh eeiaieres oe 
is applied to the cracked member as indicated in Figure 2.5, 


the steel strain at the crack is given by 


m 
" 


B/N ASE) 


pe ad as ae 
(pA, ) (nE, ) 


= P 
| p_A ) 
(1+np) 


P(1+np) 
npAgE. 


At a distance bey 12 from the crack, perfect bond is assumed 
between steel and concrete. The steel and concrete strains 


are therefore equal and given by 
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€c = P 
AeEc 


Assuming a linear variation in steel strain within the 


cracked length Ley, the average steel strain is given by 
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Therefore Ker = P/eéave 
eee 
= Kg 1 + 1 
2np 
= akg (2180) 
where a = Stiffness “reduction coefficient 


= 1 (C2847.9)) 
ae | 
2np 


Figure 2.6 from the CEB State-of-the-Art report (1982) 
shows the distribution of steel stress for a ribbed bar 
embedded in a concrete block for the working load and 
ultimate load levels. At the working load level, the stress 
distribution is close to linear and tends to become more 
linear at higher loads. The assumption of a linear steel 
Strain distribution (or equivalently, steel stress 


distribution) is thus justified. 
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2.3.2.2 The Length of Cracked Region 


The reduced stiffness of the member is affected by the 
lengtnmobethescrackedsregionibe,;toverm whichethe steelestress 
is assumed to vary linearly as described above. Several 
approaches are considered for determination of the cracked 
length Ley. It should be noted that all theories presented 
here deal with the cracking of hardened concrete and assume 
that the steel remains elastic after cracking has taken 
place (Beeby, 1979). 

1. 'No-slip' Approach 

The stress distribution in concrete within the 
cracked length can be approximated using the idea of 
Stress diffusion. It is assumed that plane sections do 
not remain plane and that bond failure does not occur at 
the time the cracks developed (Beeby, 1979). Hence, 
Phere wTSenOsSlioein the vicinityrotstncmecracks al .e. 
Lslip = 0. 

The tensile force in the steel bar at the cracked 
section is spread into the concrete by load. One may 
assume that the transfer force spreads into concrete 
roughly along a cone of slope 1/k, as shown in 
Evgurem2a,.) Them stresssoutsideathemconesiseconsicered 


ZELOMmEnrom Figure 2. a#,mtne Crackcdmenqthnars 
Ler = 2kic G2n20)) 
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face of the bar to the concrete surface. The slope of 
the diffusion cone varies within a certain range and it 
can be calibrated according to test data. The diffusion 
concept was adopted by Bazant and Oh (1983), who took 
k, = 1.4, and Clark and Spiers(1978), who took k, = 1. 
"Slip' Approach 

To illustrate the concept, consider a segment of a 
cracked member subjected to pure tension as shown in 
Figure 2.8a. For simplicity, assume that the tension 
reinforcement consists of a single bar of diameter 9. 
When the bar segment AB is considered as a free body in 
Figure 2.8b, the tensile force T carried by the steel 
alone at B must be transferred to the concrete by bond 
stress'T) over the development length Ly. The maximum 


force that can be transferred to the concrete is 
a = Aft 
If +) is the average bond stress, then 


au 


THAGLY = Auf, 
Multiply the above expression by $/4, 


TH(mo*/4)Ly = A,f+o/4 


or TpAcshlg = Ac £4 o/4 
or erp Die (2524) 
cas 


where A, = area of the steel bar 
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If the ultimate bond strength r 1s reached at the 


ult 
maximum load, then the bond stress T may be substituted 
by aes And because Tit +8 DrOopone1oOnals to f, hOnga 
given bar type (Beeby, 1979 and Leonhardt, 1977), the 


above expression can be rewritten as: 
Ly = k2¢/p OP F573) 


The above approach had been considered by Saliger 
(1936). It was assumed that plane sections remain plane 
within the concrete and that there will be slipping 
between the steel and the concrete. Bond failure is 
thus assumed to occur at each crack (Beeby, 1979). 

Since the crack affects the stresses only within a 


dustances+Uas trom thescrack,  thescrackecalengquimis 


tt 
! 


Cia 2Ly 


2k2¢/p (2.23) 


and because slipping occurs over the same distance, 


The coefficient kz can be derived from the ultimate bond 
strength suggested in the 1963 ACI Code for the 
development length. The derivation is shown in 


Appendix A. 
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3. General approach 
Beeby (1979) explained that the two forementioned 
approaches should be considered together to describe the 
actual cracking behaviour of axially reinforced tension 
members. A general formula for crack spacing is given 


by (Beeby, 1979; Ferry-Borges, 1966): 
S = kic + k2¢/p 
The cracked length can then be expressed as: 


2k cert O/B) (20'5)) 


and the slipping length is 
Lelip = 2k2¢/p : Zao) 


Note that Equations 2.20 and 2.23 are only special cases 
of Equation 2.25. Table 2.1 summarizes the different 
values of k, and kz used by others (Bazant and Oh, 1983; 
Clark and Spiers, 1978; Beeby, 1979 and Leonhardt, 


Ove e, 


2.4 Solution Algorithm 


Step ain Refering to Figure 2.2a, calculate the following 


constants for the member. 
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TOLawegsOsSsmaled, Ag 


Net concrete area, Ac 
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Steel ratio, p A</Ac 
Modular ratio, Nae Dicey fp 
Equivalent transformed concrete area, 
Ant =a ietanp) Ac 
Stiffness at uncracked section, 
Kg = AeE 


Stiffness reduction coefficient, 


Calculate the shrinkage strain at which first 
cracks forms (m=) “from Equatione2 414. 
Gheckeit "che Shrinkages strane, i Segceace: 
thane the uUltimatessnrinkagesstrai nee... 

Lio ccna an Signe emc rac hmCOES not develop. 
The analysis is completed. 

Thee 


< € the member cracks and 


shm shu’ 


member stiffness is reduced. 
Check if total cracked length is greater than 
the span length of the member. 

if mL > L, then no further cracks can 
form since the transfer length available will be 
insufficient to develop the concrete stress up 
to the concrete tensile strength and therefore 


no new cracks will form. 
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concrete stress in uncracked section after 
BOUMattOnmOLe Clack. 

Step 6. For the formation of subsequent shrinkage 
CYackSpmseteite=. 2.1.5 mm emancereneqtmste nomic) 


Om oE. 


A computer program as listed in Appendix B was written 
to perform the above analysis for given AgomAcw aC by Ley 
Ec, Eo, kan on 0) and Eshue 


In the following section results are presented of a 


Parametric study done to evaluate the effects of significant 


Parameters on shrinkage crack spacing. 


2.5 Parameter Study 


The number of shrinkage cracks that would occurs ina 
member depends primarily on the ultimate shrinkage strain, 
Span length, reinforcement ratio, reinforcing bar size and 
the tensile strength of concrete. ACI Committee 209 has 
recommended an average ultimate shrinkage strain value of 
800 millionths for concrete with 100 mm or less slump, 


and minimum thickness of members 150 mm or less, and 40% 


or less ambient relative humidity. Since the basic analysis 


procedure assumes full fixity at ends and ignores creep 
effects, an ultimate shrinkage strain of 400 millionths in 
addition to the recommended value of 800 millionths was 


considered as an approximate allowance for partial fixity 
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considered as an approximate allowance for Pantial Siva ty 
and creep. Partial fixity results from both the effects of 
flexible supports and the fact that slab systems are 
normally cast in parts (Martin, 1971). Creep reduces 
stresses induced in concrete by restraint of shrinkage. 
Thirteen members with the following basic properties 


were analysed. 


Ag = 190mm? /mm 
c = 154mm 
Ec = 25900MPa 


E, = 200000MPa 


The following parameters were varied: 
(a) Member span length, L 


(by) Steeltarean AS (mm? /mm) 


S 
(c) Reinforcing bar size, ¢ 


(i tensile strength toteconcretcmats 


To compare the results of using different approaches 
for estimating the cracked length, three analyses were 
Carried out on each member with the following variations: 


(a) Ley = Ph Odehe ky = 1.0 


(b) Lor es oe ky = 1.4 


1S) dhe 2k2¢/p, kz derived from the development 


length in Appendix A. 


Figure 2.9 shows a comparison of the cracked length 


estimated by different methods. Beeby's expression, 
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Leonhardt's expression and the development length give 
Similar cracked length and should yaeldesimivarianalyincal 
results; therefore Beeby's expression and Leonhardt's 
expression are not considered in the present study. 

Each of the parameters listed above will be considered 
in turn. A summary of members analyzed is given in 


Table 2.2. Details of parameters investigated are given in 


TaDLee Ziti 
2.5.1 Span Length 


Four members M1, M2, M3 and M4 of span lengths 2000mnm, 
4000mm, 6000mm and 8000mm respectively were considered in 
this study. For each span length, tensile stress induced in 
the slab member as a result of restraint of shrinkage is 
plotted in Figure 2.10 as a function of shrinkage strain. 
Ponsthiseseries,.bL., was takeneas 2k, copawithsk, =esl.4)q AS 
discussed previously in Section 2.3.1, a sudden drop in the 
induced force occurs at formation of a crack. The analysis 
was continued to an ultimate strain of 800 millionths. 
However it can be seen that the number of cracks formed for 
any value of ultimate shrinkage strain less than 
800 millionths can also be determined from the plot. For 
example in the case of a span length of 4000mm and ultimate 
shrinkage strain of 800 millionths results in four cracks 
forming, whereas at 400 millionths ultimate shrinkage 
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The results show that the longer the length of the 
member, the greater the number of cracks. In effect the 
crack Spacing becomes essentially the same for each span 
Yength as the ultimate shrinkage strain increases, This 
aspectwor behaviour can ibe attributedmto ther: tactethat 
formation of a crack in a long member causes less reduction 
in stiffness than in a short member and consequently a 
higher residual stress in the uncracked region as can be 


seen from the plot. 
2.5.2 Steel Area 


Members M5, M3, M6 and M7 have steel areas 0.4mm?/mn, 
0.6mm?/mm, 0.8mm?/mm and 1.0mm?/mm respectively. These 
steel areas correspond to the usual range of steel 
percentages used in concrete slab systems. Figure 2.11 
shows the results of the analyses of these four members 
indicating that the number of cracks increases with area of 
reinforcement. Again, the number of cracks formed is 
related to the effective stiffness of the member after 


iittiale cracking takes place. 
2.5.3 Reinforcing Bar Size 
The effects of bar size on the number of cracks can 


only be revealed by using the "slip' approach in estimating 


Cracked length (length of reduced section stiffness). The 
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formula for the basic development length of tension 


reinforcement as stated in ACI Building Code is 


Ld = O.019AHfY/VE"¢ 
ee 0.058¢fy 


og 300mm, 


whichever is greater, for 35M or smaller bars. The 
expression shows that the development length increases with 
increased bar diameter. 

As cracks form, the member with larger bars will have 
more reduction in stiffness because of longer cracked 
length, and thus the restraint stresses are lowered. Fewer 
cracks will therefore form in the member. This is shown in 
Figure 2.12 for members M7, M8, M9 and M10, each with the 
Same steel area but with bar sizes 10M, 15M, 20M and 25M 
respectively. The use of a larger bar with longer 
development length has the effect of lowering the induced 
concrete tensile stress and reducing the number of cracks. 
However, a longer development length implies an increase in 
Slip at the loaded end of the bar (i.e. at the face of a 
Crack hum Tooshigh "asslip may "produce=excessivercrack® width 
and thus fail to satisfy serviceability requirements for 
Grac kung’. 

The results show that for an ultimate shrinkage strain 
of 400 millionths, the use of smaller bars does not 
significantly increase the number of cracks but it has the 


advantage of limiting the crack width. In slab design, 
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using small bars such as 10M and/or 15M for reinforcement is 
good practice as far as Controllingtcrackawidth is 


concerned. 


2.5.4 Tensile Strength 


To study the effects of variations in tensile strength, 
four members M3, M11, M12 and M13 were analyzed. Member M3 
has the value of tensile strength based on the ACI 


expression for modulus of rupture, 
fr. = 0.67f¢ MPa 


Members M11, M12 and M13 have reduced tensile strength 
values of 0.50, 0.33 and 0.17V7ET respectively. The results 
of the analyses are presented in Figure 2.13. As expected, 
the number of cracks increases with decrease in tensile 


strength. 


2.5.5 Length of Reduced Member Stiffness 


The general approach and the 'slip' approach to 
estimate the cracked length L,,, which require knowledge of 
the reinforcement ratio and the bar size, is more difficult 
POMUSGMinb analysis, | Ine@cOntchaSt mticm mO=ollD Meapproccnimic 
very simple. An effort is therefore made to obtain a k, 
value in the expression L., = 2kic, such that the results 


will be in good agreement with those obtained from the 
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"slip' approach. 

Table 2.4 is a summary of the total number of cracks 
formed in each member. Each member was analysed by the 
three different procedures for estimating Ler mentioned 
eankren. §From thertable,hitgisefoundathateresultsmtorazk.c, 
ky smi. 4) eande2(Lq) jagree?welle) jAt valuenofekai =p 124)is 
therefore suggested for use in the analysis procedure. In 
most concrete slab systems, where span lengths range from 
4000mm to 8000mm, reinforcing bar sizes 10M and 15M with 
tensile strength values Om6ytcer the total number of 
shrinkage cracks formed is only one or two, as can be seen 
from the results for members M2 to M8 based on ultimate 
Shrinkage strain of 400 millionths. using the 'slip' 
approach. For a continuous slab system supported on 
columns, a single shrinkage crack would correspond toga 
crack along each column line while two cracks per span would 
correspond to cracks along the column lines and the mid-span 


lines. 


2.6 Summary 


A mathematical model to analyze the uniaxial shrinkage 
behaviour of a completely restrained reinforced concrete 
member has been developed. In the analysis, elastic 
characteristics are assumed in both concrete and steel and 
effects of creep are not explicitly included. The 


analytical procedure involves allowing shrinkage to occur 
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freely in the concrete, then applying the necessary forces 
required to maintain compatibility of deformations. These 
forces produce cracking as shrinkage strains increase. The 
member stiffness is reduced when cracks form, such that the 
internal force decreases by some extent and increases again 
as shrinkage continues. The reduction in stiffness is 
assumed over a cracked length within the vicinity of 
influence of a crack. The reduced stiffness in the cracked 
region is developed by assuming a linear steel strain 
distribution over the length concerned. Three approaches, 
"no-slip' approach, 'slip' approach and a general approach, 
for estimating the cracked length are described. 

A parametric study was performed to investigate the 
effects of several parameters on the total number of cracks 
forming in members subjected only to uniaxial shrinkage. 

The study illustrated the effects of variations in span 
length, reinforcing area, reinforcing bar size and tensile 
Strength. The results show that the number of cracks formed 
is sensitive to the variations in concrete tensile strength 
but relatively insensitive to variations in other 
parameters. The results also indicate that for span lengths 
of 4000mm to 8000mm, one or two cracks per span length can 
be expected under normal circumstances. A value of k, = 1.4 
used in the 'no-slip' approach was found to give results 


Similar to the 'slip' or general approaches. 
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Development Length OM) Cie Omnia 0.627 F2* 


0.12 0.50/F 
0.082 0.337FE 
0.041 Desi eS 


Ape = 43 0MPa 


Table 2.1 Summary of values of ky and k» 
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Member 


All members 


Tabpler2.c 


Steel area Bar size 


(mm? /mm) 


A. = 190mm2/mm 


g 

6 25900MPa 
ec = 154mm 
a 8S Poe 


Summary of members analyzed 
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As = 0.6 mm2 /mm 
10M Bar size 
fi = 5.415 Mea 


Span length M1 
M2 
M3 
M4 
Steel area M5 
M3 
M6 
M7 
Bar size M7 
M8 
M9 


6000 mm 
10M Bar size 
fr = 3.415 MPa 


L = 6000 mm 
Y.0 mm2 /mm 
34105 MIRE 


> 
10) 
Il 


Tensile strength 3.415 MPa 6000 mm 
ees ThEY 10M Bar size 
183 22 As = 0.6 mm*/mm 


Oro tal! 


Table 2.3 Summary of member properties in each parameter study 
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Shrinkage strain versus time curve 
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(a) Fixed support member 


(b) Release one support 
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(c) Compatible member 
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Figure 2.2 Symmetrically reinforced slab element 
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Figure 2.3 Member response curve 
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Reduced axial stiffness 


(a) Fixed support member with one crack 


(b) Release one support 
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Figure 2.4 Cracked slab element 
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Figure 2.5 Strain distribution in a cracked member 
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Figure 2.6 Distribution of steel stress along 
the anchorage length 
(Cis SIs) 
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(a) Slab element subjected to pure tension 


(b) Free body diagram of bar segment AB 


Figure 2,8 ‘Slip’ Approach 
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3. ANALYSIS OF SLAB SYSTEMS INCLUDING EFFECTS OF CRACKING 


3.1 Introduction 


In the analysis of reinforced concrete slabs, effects 
of cracking due to both shrinkage restraint and transverse 
loads should be considered. Cracking reduces the overall 
flexural stiffness of slabs, resulting in increase in 
deflection and redistribution of moments. Serviceability 
may be adversely affected by excessive deflections and 
extensive cracking. 

The finite element method provides a means to take 
account of different properties in each region of the slab, 
in particular, the reduced flexural stiffness due to 
cracking. In the finite element model described 
Subsequently, cracking due to transverse loads and due to 
restraint of shrinkage is modelled. Long-term deflections 
due to shrinkage curvature and creep are also included in 
the model. 

Analyses are carried out on three reinforced concrete 
slabs to study the effects of number of shrinkage cracks and 
variation in modulus of rupture. Results of the study are 
used to develop a simplified method to calculate slab 


deflections including shrinkage effects. 
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3.2 Finite Element Model 


3.2.1 Assumptions and Limitations 


The finite element computer program SAPIV (Bathe, 


Wilson and Peterson, 1974) was modified to account for 


reduced stiffness due to cracking, using a procedure 


proposed by Scanlon and Murray (1982). This procedure 


contains the following assumptions and limitations: 


ile 


Linear elastic response is assumed for reinforcement and 
for concrete in compression. The analysis is therefore 
not valid for the post-yield range of behaviour. 

Effects of cracking are included using the effective 
moment of inertia procedure proposed by Branson (1963). 
This procedure provides for a transition in flexural 
stiffness between the uncracked and fully cracked 
limits. Originally proposed for beam analysis, the 
procedure is geherathizeaiein this investigation to 
consider two-way (plate) action. 

The tensile strength of concrete under biaxial stress is 
assumed to be equal to the uniaxial value of modulus of 
rupture. Figure 3.1 shows the biaxial failure envelope 
of concrete in terms of principal stresses. Combined 
tension and compression loadings reduce the tensile 
strength; however, in the service load range, the 
reduction is not significant and therefore the use of 


uniaxial value is reasonable in two-way slabs (Scanlon 
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and Murray, 1982). 

4, It iS assumed that cracking is initiated when My or My 
exceeds the cracking moments Mcy. My and My are defined 
for the local coordinates of the plate bending elements 
aSmSnOwn glngnigure 3.2) tinegeneral, stor thes finite 
element mesh layouts eee in this Study, pthe tocar 
coordinates coincide with the slab global coordinates 
defined by orthogonal reinforcement directions. 

5. The shear modulus, .G; efor concrete; required for 
determination of the plate torsional rigidity is assumed 
to be unaffected by cracking. In another study (Hand, 
Pecknold and Schnobrich, 1973) it has been shown that 
plate bending analysis is not significantly affected by 


variations in the assumed value of G. 
3.2.2 Flexural Stiffness of Cracked Region 


The analysis of a reinforced concrete slab is treated 
aSma oroblem of orthotropic plate bending.= the plane stress 
constitutive relations for orthotropic material are given 
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The moments are related to the curvatures by 


Mx Beane vB he e p x 
T2Ui-vyry) TaCi-vyv,) 

My = Vv Bene By nh? e Py 
(TACSTar a Gee 

Mx y ® e ae Dim, 


re) 


My and My can be obtained for each element from an 
analysis of the uncracked slab and then checked against Moy 
Bom cuhacking. If cracking 1s detected, reduction) in 
flexural stiffness in each direction can be accounted for 


using the Branson's expression (1963): 


aca (Mop /My) *Tg +o ia Mer Me aa eles (e303) 
Tey = (Mop /My ) *Tg 1 [ a (Mcp /My ) * ] Tory (3.2) 
where Moy = ErIg/Yt eles) 


The reduction in flexural stiffness due to cracking is 
implemented by modifying the plane stress constitutive 


properties as follows: 


By = ayEe ve tae etal 
By = ayEe Vy = ayv (316) 
where, ay = Ig,/Ig and ay = Tey/T (S78) 
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The analysis is repeated using the reduced values until 
the results converge. Three iterations appear to be 


adequate for the slab systems considered in this study. 


3.2.3 Finite Element Model of Shrinkage Cracks 


To study the influence on load-deflection response of 
cracking due to shrinkage restraint, concrete slabs were 
assumed to be precracked as a result of shrinkage restraint, 
prior to applying transverse loads. The following procedure 
was used: 

1. Reduced flexural stiffness of precracked region 

In the analysis of the symmetrically reinforced 
member subjected to shrinkage, described in Chapter 2, a 
linear variation of steel strain within the cracked 
region Lcy was assumed. For analysis of slabs in 
bending, a linear variation of curvature within the 
precracked region was assumed in deriving the reduced 
flexural stiffness. 

At the gross section 
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and at the cracked section 
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transformed section. It is implicitly assumed that the 
ermack width is initially zero for precracked region so 
that an application of bending moment compressive 
Stresses can immediately be transferred across the 
crack. 

The average curvature within the precracked region 
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is 


ave = ( og + bcy)/2 
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If Inx is the reduced moment of inertia of the 


precracked region in x-direction, then 
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As for the cracking model, material properties in 
the constitutive relation are modified for precracked 


elements as follows: 


aie 7 Gryeats Vx = px? 


By = Aap yEc Vy = ap yY (237-51,2,) 


The analysis then proceeds as for an uncracked 
slab. 
Location and width of precracked region 

It is observed that elastic flexural tensile 
stresses due to transverse loading are highest along the 
interior column line and then along the mid-panel as 
indicated in Figure 3.3. When tensile stresses due to 
Shrinkage restraint are superimposed with flexural 
tensile stresses, cracks are expected to form first 
along the interior column line when the resultant 
stresses exceed the tensile strength of concrete. When 
combined tensile stresses are high enough, cracks will 
form next along the mid-panel. 

In the parametric study of the uniaxial member, 
the width of a cracked region of reduced stiffness was 


determined to be approximately 


=m 2 (ia) c ones) 


This relationship was used to specify the width of 


zones of reduced stiffness due to pre-cracking and 
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formed the basis for specifying the width of precracked 


elements. 


3.2.4 Long-term Effects due to Shrinkage Curvature and Creep 


Deflections of reinforced concrete slabs increase with 
time. The additional deflections are caused by creep and 
shrinkage curvature. This additional inelastic deflection 
increases at a decreasing rate during the time loading. 
Procedures recommended by ACI Committee 209 (1971) and 
435 (1966) were used to calculate long-term deflections in 
addition to immediate deflections. 

1. Shrinkage curvature 

Concrete shrinkage in unsymmetrically reinforced 
members, aS in the case of slabs, causes a nonuniform 
Strain distribution and results in curvature of the 
member. In flexural members, shrinkage is more 
restrained at the tension face where heavier 
reinforcement is placed. Therefore, shrinkage 
curvatures will have the same sign as the curvatures due 
to transverse loads and consequently increase the 
deflections. 

Expressions for shrinkage curvature are provided by 


Branson (Branson, 1963; ACI Committee 435, 1966): 
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for (p-p') < 3%, which is generally the case for 
reinforced concrete slab systems. Note that p and p' 
are in percent. 

In the modified SAPIV program, shrinkage curvatures 
were treated as equivalent temperature curvatures for 
each plate element. Deflections were then computed 
based on the equivalent temperature curvature. 

Creep 

ACI Committee 435 (1966) recommended a procedure 

for computing the deflection due to creep based on the 


work of Branson as given by 


Acp = kpCy (Aq)p 
where kyeee=) CompressilonMsteelmbactor 
~ OFB5 
150 py 
Cree CECeD coefficient 
and (Aj )p = instantaneuous deflection due to all 


sustained load 


This is equivalent to using a modified modulus of 


elasticity in the computat ionsofselasticedet lection: 
Ech e= Ec 


In the finite element analysis, the values of k,;y in 
x and y directions are obtained to compute the modified 


modulus of elasticity in the corresponding direction. 
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The reduced values are used in the constitutive 
relations to calculate deflections. In this study, 
Ky = 0.85 was used, since only small amounts of 


reinforcement are provided in compressed concrete zones. 


3.3 Finite Element Analysis and Parameter Study 


3.3.1 Slab descriptions 


To study the effects of precracking due to shrinkage 


restraint analyses were carried out for no precracking, 


precracking along column lines only, and for precracking 


along both column lines and mid-panel lines. The slab for 


which no precracking was assumed, was analyzed for a range 


of modulus of rupture values. Details of the slabs analyzed 


are as follows: 


ila 


A square two-way Slab S1, designed by the direct design 
method, was used to study the effects of variations in 
modulus of rupture on the deflection of a slab with no 
pre-cracking. The slab has three 6m span in each 
direction, a thickness of 190mm and 550xX550mm square 
columns. The slab is a flat plate with no edge beams or 
drop panels. The thickness selected satisfies minimum 
thickness requirements of the ACI Code. Column size and 
reinforcement were selected on the basis of strength 
requirements. The structure layout is shown in 


Figure 3.4 and a finite element mesh layout is shown in 
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Pigureme.5)% ‘Because of symmetry, only a quarter of the 
Slab need be considered in the finite element analysis. 
Results of the analysis are presented in Figure 3.6 in 
the form of load-deflection curves for the center of the 
exterior panel. The deflection at this location was 
observed to be the maximum for all values of modulus of 
PUpEUGEe.. Wehrgurel3 /ieéshows tthe) addita onal adem@lections 
due to creep and shrinkage curvature. 

The same slab system appears in Figure 3.8 as slab S2, 
except that the slab is precracked along the interior 
column line. The distance from extreme compression 
fiber to the face of tension reinforcement c is assumed 


to be 154mm for all sizes of reinforcing bars in both 


directions. The resulting width of precracked region is 


2(1.4)154 


i 
iT) 


Gir 


R 


440mm 


Figure 3.9 shows load-deflection curves for the 
center of the exterior panel. 
Slab S3 1s the same slab system but precracked both 
along the interior column lines and mid-panel, as shown 
wimbedures3. 10... Thiseslabsrepresentsaunesmoremsevere 
case where tensile stresses induced by restraint to 
shrinkage are high enough to produce cracks at both 
locations. Figure 3.11 shows the load deflection curves 


for center of exterior panel. 
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3.3.2 Parameter Study 


Bttectssoseprecracking 

Moments along the interior column strip of 
precracked slabs S2 and $3 are compared to that of 
uUnerackedsstabeSi, in Figuremscmaies mS labeS2anas the 
highest moment at the mid-span, which is a result of 
redistribution of moments when a shrinkage crack is 
formed along the column lines. The cracks at mid-panel 
of slab S3 reduce some mid-span moments and redistribute 
the moments to the columns. The significant effect of 
precracking on the redistribution of moments is 
illustrated here. 

The effect of various degrees of precracking on 
deflections is shown in Figure 3.13. The excessive slab 
deflections of slab S3 would probably occur under 
extreme adverse condition and therefore represents a 
conservative deflection calculation. Slab S2 represents 
an average condition where the deflections do not 
deviate greatly from the uncracked slab $1. 

BEfects of variation in modulus of srupture 

For slab $1, deflections were computed based on the 
effective moment of inertia using different values of 
modulus of rupture in Equation 3.5. The modulus of 
rupture f, in Equation 3.5 may be replaced by the 


effective modulus of rupture f,: 


Mer = felg/Yt (3.14) 
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The effective moment of inertia Ie is then 
calculated using the reduced value of Moy and subsequent 
deflection calculations are carried out in this 
Simplified method to take account of restraint stresses 
(Scanlon and Murray, 1982). 

Variation of deflection with effective modulus of 
rupture at different load levels is shown in 
Figure 3.14. It is observed that at lower load level, 
variations of modulus of rupture has little effect on 
Slab deflections. Indicated in the plot is also the 
deflections of the two precracked slabs. The results 
Show that effective modulus of rupture values of 0.3 to 
0.4¥£¢ MPa at service load level corresponds 
approximately to the case of precracking along column 
lines only. Effective modulus of rupture of 0.2 to 
Oe sy ie applied to an uncracked slab corresponds to the 
case of a slab precracked along both column and 


mid-panel lines. 


3.4 Summary 


This chapter describes a finite element analysis of 
slab systems that may be precracked in prescribed locations 
as a result of restraint of shrinkage. The analysis was 


applied to a slab system under no precracking, precracking 
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along column lines only and precracking along both column 
and mid-panel lines. These cases represent the range of 
conditions likely to occur under none EOndi tt ons pany, 
means of a brief parameter study it was shown that the 
effects of precracking due to shrinkage restraint could be 
accounted for by applying a reduced effective modulus of 


rupture in the analysis of an Bnitially uneracked slab. 
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Working 
Load Range 


0.5 
cou 


Failure 
Envelope 


Tension Positive 


oF 0, ore principal stresses 


Figure 3.1. Biaxial strength of plain concrete 


(Kupfer, Hilsdorf and Rusch, 1969) 
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Figure 3.2 


Plate bending element 
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High flexural stresses along 
interior column lines 
and mid-panel lines 


Continuous edge 


Figure 3.3 Tensile stress distribution in slab along 
column lines and panel centre lines due to 


uniform transverse load 
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All columns 550mmx550mm 


Figure 3.4 Structure layout of slab system 


Figure 3.5 Mesh layout for Slab Sl 
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Figure 3.6 Load-deflection curves at center of 
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LEGEND 
aa Immediate deflection 
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Figure 3.8 Mesh layout for Slab S82 
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Figure 3,10 Mesh layout for Slab 53 
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4. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


4.1 Summary 


In continuous reinforced concrete slab systems, the 
effects of shrinkage on slab deflections are often 
underestimated. The current ACI Building Code considers the 
additional deflections due to shrinkage as a result of 
shrinkage warping which arises from nonuniform shrinkage 
through the section. Reinforcement of different amounts in 
the two faces of a member, which is usually the case for 
reinforced concrete slabs, is the main cause of nonuniform 
shrinkage. The Code however does not explicitely recognize 
the induced tensile stresses in concrete due to restraint of 
shrinkage in continuous slab systems. These induced tensile 
stresses may contribute to cracking of concrete resulting in 
an increase in deflections. This study presents a procedure 
using the finite element method to estimate the effects of 
restrained shrinkage on two-way slab systems. 

A new mathematical model is developed in Chapter 2 to 
analyze the uniaxial shrinkage behaviour of a completely 
restrained reinforced concrete member. Cracking occurs when 
the induced tensile stresses in concrete exceed the concrete 
tensile strength. Reduction in member axial stiffness is 
made by assuming a linear steel stress distribution within 
the zone of influence of cracking, or cracked region. The 


length of the cracked region is based on a crack spacing 
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formula. The analysis determines the number of cracks for a 
given length of member and the shrinkage strain at which 
each crack forms. A parameter study was carried out to 
studvetnie ertectsmotevariations inespan length, steel area, 
reinforcing bar size and tensile strength on the total 
number of shrinkage cracks forming. Results obtained are 
used to. select the locations of precracked regions in the 
analysis of a slab system subjected to restrained shrinkage 
and transverse loads. 

Chapter 3 describes the finite element analysis of a 
two-way reinforced concrete slab system subjected to 
restrained shrinkage and transverse loads. The effects of 
restrained shrinkage are modelled by forming shrinkage 
cracks prior to analyzing the slab for applied loading. 
Precracking is accomplished by reducing the flexural 
stiffness of elements at the locations of shrinkage cracks 
based on a linear curvature variation across the width of 
the precracked elements. The effects of the number of 
shrinkage cracks and variations in the modulus of rupture 
are studied. Results of the study are used to assess a 
simplified deflection calculation procedure proposed by 
Scanlon and Murray (1982) for slabs with restrained: 
shrinkage. Scanlon and Murray proposed values for reduced 
effective modulus of rupture based largely on intuitive 
grounds. The present study has provided a basis for 
specifying values for effective modulus of rupture based on 


numerical studies including explicit consideration of the 
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presence of shrinkage restraint cracks. 


4.2 Conclusions 


The following conclusions are drawn from this 


investigation: 


1. Uniaxial shrinkage analysis 


a. 


Reduced axial stiffness in cracked region is 


expressed as 


where re EI 1 


The 'no-slip' approach can be used to approximate 


the length of the cracked region: 


Ler = 2kac 


Using a value of k; = 1.4 in the above expression 
yields similar results as those using the more 

complicated 'slip' or general approaches. Having 
the advantage of simplicity, the stress diffusion 
expression with value of k, = 1.4 was also used in 


the two-way slab analysis. 


2. Reinforced concrete two-way slab system 


a. 


Shrinkage cracks are modelled by precracking 
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elements at the locations of the cracks. The width 


of the precracked elements is 
Lop = 2(1.4)c 
and the reduced flexural stiffness is 


(EI) cp = apEclg 


where 


i) 
N 


9 


b. Based on results of a brief parameter study, a 
reduced effective modulus of rupture in the range of 
On 2eto 004) i scanibesusedstosacccuntmrormertiectcmor 


restraint stresses on slab deflections. 


4.3 Recommendations for Further Research 


The uniaxial shrinkage analysis can be extended to 
include the following features: 
1. effects of creep (see Appendix D), 
a unsymmetrical reinforcement placement, 
3. varying amounts of steel area along the span, 
4. presence of flexural stresses, and 
5. supports with partial fixity (see Appendix C). 
A complete parameter study for the extended analysis 


can then be carried out. 
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The finite element analysis of two-way slab systems may 
be extended to include the following capabilities: 
1. Consideration of post-yield behaviour; 
ene Modi ficationgofethe.crackingscriterionstosconsides 
principal moments that may be at a different orientation 
than the element coordinate system. 
Further research may be directed to a detailed 
Parameter study of factors such as the effects of drop 


panels, edge beams and depth to span ratio. 
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APPENDIX A 


Coefficient kz for Development Length 
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ayy 


Brome Equations #2221 and 2.22, thescoetticientak, =tor 


the development length Ly is given by: 


K2 = fr 


mn 
4 Tuylt 


for pure tension. 


The average ultimate unit bond stress capacity ry} is 


given in the 1963 ACI Code as: 


71 tee ana < 800 psi 


or Z0Y eae =O OMe 
? 
for #11 and smaller bars, which are usually used in concrete 
slab design. Using this value of ry, ]4, kz can be computed 
for different values of modulus of rupture f, as summarized 
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f, (Mpa) 


Summary of values of ky for development length 
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Computer Program for Uniaxial Shrinkage Analysis 


89 


Ps : = : / >be : > iF ear . 
eke) Reid © se: écr 7 
Vegas Meee oe me Se 


a ey : - 
ne ee ae 
7 _ | 


90 


* 3} 
YyOHuae, paxyoseud [e2}0O} BY} YOSUD - pv d3IIS * @) 
* e) 
eK KR OK OK 2) 
2) 
666 OL0D (NHS3°1D°HS4)4I 
3) 
ORK KK OK OK OK 2) 
* oO 
uUpeuzyS SeHeyutuys ey HY98UD - E€ dJIS * 3) 
* 2) 
OK KK KOK OK 0) 
3) 
(OILVY * GQIOH) /y4*43aV = HS4 OL 
SA = daioan 
OnE =) OFLLVa 
bt = W 
re 
KKK ROK KO 9) 
* 2) 
uULeu}S eHexAuULUYS BY} BeELNOLeD - © d3Z1iS * 3) 
* 3) 
2K KOK OK KOK k 2) 
&) 
(tO07‘S9)3SLIYM 
NHS3A‘' SH EM VII‘'S3A ‘OHS SV‘'44‘°OR OV VO DV ATILIL (0007 ‘9)3LTaM 
2) 
diiS7 + O03 Ct D4tw = y4o71 
O3°@x OHY/VIG*«tH = dI1S1 
DH «x ((OHY*Nx?)/OIJ0°4++O090°4+)/O30°1 = Yon 
eke = oP | 
OV* (OHY*xN + OFO'L) = dav 
03/S4 = N 
OV/SV = OHd 
SV-DV e= OV 
2) 
ORK KK KOK OK 3) 
* 2) 
“deaqQueaewW Stj}, UOJ S}PUCPSUOD SY} SPELNIOLSD - | d3lS * 9) 
* 2 
KK KK OK OK 3) 
3) 
ROR RK OK KOK OK KK stsA|euy A KK KK KO 
9) 
NHS3‘VIG' cH‘ tH‘ TO * 
“a4‘'Sa‘'Oa'SV'DV * 
“A1LIL(L66=G0N3 ‘OOO! ‘G)avaa Ss 
3) 
o 
eyeqg yndul xx**xd 
2, 
44dILS IVNYSLXS 
Yon’ DH'YDT' 1 /dOY¥d/NOWWOD 
(O@)ATLIL bxdFDILNI 
G10 * MANX‘ YOM OH'N' dI IST YO CM EMS TI WV3a 
ab 
2 
-uequew padu0oyULeU 6) 
ALLeOtuzZewwAs e 4oOo asucodsau esHe7nHuruys ey. saZA,~eue weuboud siys 2 
5) 


“VOtNONAD 


ittat * cnn i . iy Dae a 
* hea 
: oe 7 7 7 7 
| =) Ou 
iMeren, CM A1G) soa ase: on 5 Lee 


_ 


whee! semi vie ar slelee ne 


= ¢ aa i? RK a a - 
: _ 
| 


(CLAW © DI0MY HT 


_ 8: 
= * 


.a 
— 
= 2 
sf 


= 7 
ube» ape dri ote 4 a a 
a ece SGD TSS re font! 


3 aeons 
: ° 
7 
‘tee! neice: (eshs et) Gh ~ & an be 

: * 


? 
: 
7 
. 
7 - 
7 


sal 


///v. Glia v SIDS b Sha 7 ; 


* 
/,Ae}owetq ; * 
APO Sete O1}eu [eas eeuy 18ers Fig Ui * 
/;VIVLAG LNAWSOYNOSNISY, ‘, , * 
Ff Ves D Ge aS eee * 
/,4tS OLtsuay Aytor1zyselg Badly eyeJ9U0D 7-7) ; * 
raislexe| elitereyNebeiyn ap , * 
Wi he GLNS? 5 7 x 
/,4}6ua7 ueds 4BAOD xXeW eauy Sssouy pte er 
/,11VL50 S,,/YSEWAW .///vVOC', »/',b.)LVNYOS OOOT 
(O° O149/0° OL 4G/vVO7C)LVWYHOS OOO 
KK ROK OK KOK K HOO LVWYOS 2K KK KK kk 
dO!lS L466 
G O105 
(OLOE*9)3LIYM 
(OCOE*S9)ASLIYM SEE 
GS 0109 
(OLOE‘S9)3ALIYM 
(OOOE*9)3S1TYM 666 
2K KOK ROKK KK NOITLVNIWYSL KK KK OK KK 
Ot OLOD 
MANH = GQ10H 
+ +W=W 
2K OK OK KK OK OK OK 
Sd 
syoeud JUaNbsesqns. uvOos SuNnpaedoud ay} }FeEaeday - 9 daIS x 
* 
EK OK KOK OK OK OK 
yO1L°YOS4A‘9DS4°O4*HSA‘W(OTO7T* 9) ALTIUM 
4ndino 
OHY/94 - YOSd = 9S4d 
SV/1L4Nd = awOS4 
OLTLVY+HS3*S3%*SV - d = L3ANd 
Sessau}S |8a}S Bay azyepNoLe|) 
avV/d = we 
OILVa * HSS *« MINH = d 
V/(dI1S1*W-1) = OILva 
(W)44ILS = MIN 
3K KK OK Ok Ok ok 
* 
SessedysS SY} SPELNIL2D - G dwJI1S «x 
* 
2K KK KOK KOK OK 


S66 O109 (7 1D ADIL) SAI 
YO1+W = AOTL 


De OK te oe ae Ook oe 


0) 
3) 
8) 


0) 
0) 
0) 


OlOlOnO © CIC OLOIS 


OOO 


COOUO0U00 


90 


‘beeenss ‘teets ——- 10 


orraneizae "eae 2s 
Za\irw 
SND = RA ~ 
7 7 7 : ae 
TS. 
nt AOE ET 


a - — 


PE ie) ae piwelapegah om wo oH 


7 Pi ae 


ol spaapsuyss “aie smears 


vi ances ‘ITs renee : c eevee 


nae — ou -_ 7 
ek a pint 


. SATIN 2°. 2a8 uIP 1 Serate... 7 ead cdl y , 
ingne) een? jevou “ah Per =< ee 
; (ia: Be _ 
‘) JPeTSO : 7 29y¥p2" oe 
they yoiaiaaat3 ir tee ya . i ; 
Le 4 bea i 7 
vi URNS TaREanealeaa a 


ad at og 
{i S\*on4 Wee's eB es 4a°9 reare an - 
“ 


wWTee eis 


92 


erty sO Pug 


GN3 St 

Naniad OSt 

(YON/YOT4W + DH/(CYOT*W-1))/7 = 4SILS 6vl 

2) BV 

Yon’ DHS YO1' 1 /dOYd/NOWWOD Ly 

Yon DH’ YNOT' 1 Wad Ort 

3) Gtr 

AK OK KK KOK OK KOK OK KOK OK 0) prt 

* 2) Evt 

“SyoeUD YW YIM SSOUSF}RS * 5} cv 

eGeusrae s,uequew ey} Sse}yen,eASe uUol YOUN; SUL os 9 ‘wea 

* 2) Ort 

SOKO OK OK RK ROK KOK ok 3) 6El 

2) BEL 

(W)44ILS NOILONN4S WV34 Let 

@ SEL 

2K EKKO ROK OE RRO RRO ROR ROKR I I OK aOR EI IE IK OK aK aR KK aK a aK aK ak aK aK ak ok ak ak GE 
ON3 vel 

@) ES 

(,**«* Y}GHUe, uvequeaw speasoxe yyOHua, paxYoeud [e301 xx ,/)LVWYOS OZOE cel 
(,*x** W4OF LLEM S4HDCUD GAYUPUNY ON xxx ,/)LVWYOS OLOE bet 

(,**k* SNLVA BPeWIYLN spssoxs UlLeuYS BHeyuLYYS xxx ,/)LVWYOI OOOE OEt 
(ct Sta’ vp GhDOE* XG’ bh GLA‘ LI)LVWYO4 O7OT 6c 

(/,4;6uey ~ 40 ; x Bet 

*, (49249) 98a31S (ssoub)1eais SS38d1S DUO) , ~  ® ETAL ' 
eo abe yut sys WOBbY ga 4 * 9Tt 

PE ES WS A OW WON AY ey) Wy re a) IN de eel dal SS al AA ae ee TE ISP lo oe Gch 
JLVWHOS 410072 vo 

(///¢ Gla, = ULeuwzZs eGenxuLaYys Sez,eEWwIZIM.*, -; * Ett 

SAL UPA Hoh Fees aaah = Ne ewe at So | yap aiey. * Car 


f-° ABOO HLONSAT GaNHoOvAdD.*., , * Lee 


a ———eeeeeeeeEsEeEeEOEOOrO 


Y 
oe ae 


7 ry 


> a a 

F a v Raith 

nm - “sreret seamed eee 
d - 


CLenn RA Rea Eee rene ee Peeetreeterrerrritti rs aves slats 


me 
7 
a “4 


nenaueie 


a a 
SyOve 2 serwen wr ahiepten me pet aeiee eis 


eee 
* fener o oe 
— 


. oo oh Sie se 


APPENDIX C 


Partial’ Fixity 


a3 


94 


A simple case of partial fixity may be included in the 
uniaxial shrinkage analysis as follows. 

Consider the simplified slab system with partial fixity 
in Figure C.la. The slab is fully fixed at one end and 
attached to a column at the other end. The slab and the 
column have stiffnesses K, and Kom nespecoively:. 

If the restraint at B is released as shown in 


Figure C.1b, then the slab shrinks freely by an amount 


A = e pL CCcaly 
where e€, iS given by Equation 2.16. 
When the restraint is reapplied, a tensile force P is 


induced in the slab and the slab elongates by an amount 


as shown in Figure C.1Ic. 

The slab thus actually shortens by (A - Acjap). The 
column also deflects by the same amount as shown in 
PergquremG.c. 


Therefore, 


P SK Ameo ne Ge (C3) 
Substitute Equations C.1 and C.2 Intoeces 
P = Keleml - PL/Km)/Lc¢ 
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Keys Stiffness of the system 


L 
pairs ay an 


Ky 1S given by Equation 2.4 and Ko 2S given in 
Figure C.3 for several cases of support conditions. 


Substituting Keys for all occurrences otek sine thems prewlous 


equations in Chapter 2 will include the effects of partial 


fixity in the uniaxial shrinkage analysis. 
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Column stiffness = K. 


(a) Slab system 
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(c) Restraint reapplied 


Figure C.l Slab system with pyeyeep artless <allinny 
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Figure C.2 Free body diagram of the slab and the column 
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(c) One end fixed, one end hinged 


Figure C.3 Column stiffness 


98 


APPENDIX D 


Alternate Solution Algorithm including Creep 


che 


ae 


_ 


se a 7 
o6s2 ontbelso: witivepia nebwioe® stanse2té. 


oe 


100 


A time-incremental analytical procedure may be used to 
approximate the effects of creep in the uniaxial shrinkage 
analysis as follows. 

The quantities which must be known for the analysis 
include: Agree Acgwmacvely ty, Ec ik ye komonde OG. gel teiseascumed 
that 
Wee shrinkagerand creep start at time: t=0, 

2. the concrete has attained its design strength before 
creeping occurs, 
3. the slab element is initially uncracked, i.e. m=0, and 


4, the concrete is unstressed before shrinkage begins. 


The following constants for the member shown in 


Figure 2.2a are evaluated initially: 


Neteconcreretareay Ac = Stigie.Ac 

Steel ratio, p = A,./A; 

Cracked length, ice = 2( kee Kooy 0) 
Slipping length, Lelip = 2k2¢/p 


During any time step t, the following procedure is 


used: 

See Deen. Compute the modified modulus of elasticity due 
to creep, eho and the unrestrained shrinkage 
strain, €.sp, at time t using the formulae 
suggested by ACI Committee 209. 

Ste ue2. Compute the following time-dependent quantities: 


Modula meracl Oye = E./Ect 
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Equivalent transformed area, 
Ag = (i+np)A, 
Stiffness at uncracked Section, 


Kg = Ae Ect 


Stiffness reduction coefficient, 


a = 1 
lesa 5 | 
2np 


Stiffness at cracked section, 


Kor = akg 


Member's average stiffness, 


K = L 


mM 
(L-mbe,)/Kg + mbcgy/Key 


Member's average shrinkage strain, 


Em = é€sp (L-mL,14))/L 
Seeors. Compute the concrete stress, 
fom = Km &m/Ae 
Step 4. (a) If £..< £,, the member does not form a new 


crack. Repeat Steps 1 to 3 using the same m 
value butva new time-step, ttAt. 

(b) If f£,,2 £,, then the member cracks. The 
time at which the member cracks may be 
obtained by using the bisection method. 

The required initial bracket is 
ftiboact toelewhere ot Jo gistthe previous 


time-stepvand t"o els thescurrent jrimesstep 
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as shown in Figure D.1. An iteration of the 
bisection method consists of three steps: 


1) compute the mid-point 


2) evaluate the concrete stress f,(t,) at 
time t, by repeating Steps 1 to 3. 

3) “determine ‘a new antervareit (eee te oe 
according to the rule 

se sey ie eink fae ion Beau 


Cee ei a = 


Lope  eOtnerwuse 


The iteration continues until one of 


the stopping criteria is satisfied: 


1) | te =. Cy | Sa Or 
22) ee Peis ai Ol eet OL 
where TTOL = The acceptable tolerance in 


time, 
FTOL = The acceptable tolerance in 


concrete stress. 


After) thestime at) crackingerr ssececenmined in 
Step 4(b), increment the number of cracks, m, 
and repeat Steps 1 to 4 at time t,;, for the 


residual concrete stress after cracking (point B 
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in Figure D.1) and at subsequent time-steps. 
Step. 6. The analysis should be continued until the 


ultimate specified time limit is reached. 


A computer program which includes the effects of creep 
and partial fixity is written in Appendix E. 

The thirteen members in the parameter study is 
reanalyzed to examine the effects of creep. The modified 
modulus of elasticity and the unrestrained shrinkage strain 


in the study are determined as follows: 


Ect = E 
Cine 0 38CH) 
peice tenes 
LO eee 


€Esh = | t | 740x107 
BD mcaeet 


where Ct 


Peorecndmc.peane plotted in Figures, D- 2¢andep.seaceamfunc tion 
of time. Results are presented in Figures D.4 to D.7 for f, 
versus log t, also in Figures D.8 to Dvlietor 1 gversushe.. 
Table D.1 summarizes the total number of cracks formed in 
each member. The results show that shrinkage cracking 
occurs at an early age and the total number of cracks formed 
is close to that obtained from the previous analysis using 


an ultimate shrinkage strain of 400 millionths. 
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Concrete stress versus time curve 
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Figure D.2 Modulus of Elasticity versus log time curve 


106 


Shrinkage Strain 
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Figure D.3 


2 5 2 5 2 
ee WOE ee 


Time, days 


Shrinkage strain versus log time curve 
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Figure D.4 
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Concrete stress versus log time curve for 


errects Of Varidtions in®span length 
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Conc Stress, MPa 
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Concrete stress versus log time curve for 


effects of variations in steel area 
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Figure D.6 
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Concrete stress versus log time curve for 


effects of variations in tensile strengcn 
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Figure D.8 Effects of variations in span length 
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Figure D.9 
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Figure D.10 Effects of variations in bar size 
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Principal Moment Analysis 
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In the finite element analysis in Chapter 3, cracking 
was assumed to occur when M, or My exceeds the cracking 
moment Mcy. However, experiments have indicated that cracks 
are initiated normal to the major principal moments 
(Lenschow and Sozen, 1966). There may be elements where the 
major principal moment exceeds Mcy but not My or My. In 
such cases, cracking in the element is undetected and may 
result in an underestimation of deflection. 

The slab system in Chapter 3 is therefore reanalyzed to 
check for the principal moment in each elements. The 
principal moment is computed using the Mohr's circle as 
shown in Figure F.1. Figure F.2 shows the major principal 
moments in Slab S1 after three iterations of modified SAPIV 
analysis. Because of symmetry, the principal moments are 
shown for only half of the elements. 

The results indicate that only four of the elements has 
undetected cracking; also, all the cracked elements have My, 
or My close to the major principal moment, presumably 
because of relatively low twisting moments My y along the 
column lines. The deflection is therefore not significantly 
affected by the assumption of checking M, and My COL 


cracking in the slab system analyzed. 
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Twisting 
* Moment 


(Mo vee) 


MO M1 Moment 


(My. My) 
(M,+M_) 9 5 
Diameter, d = se = a + My 
(Mx+My ) 
M, = rs d 
(M, +My) 
Mo = Te. 


Figure F.1 Determination of principal moment 


using Mohr's circle 


All units in kNmm/mm 


MS 20.5 kNmm/mm 


* Undetected cracking 


Figure Fo2 Major principal moments in Sikeley AM 
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APPENDIX G 


Conversion Factors 
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Imperial — SI 

ee uner=—e25).-2 ei 

ie eee = 645701 G2 mm 

1 lb(mass) = 0.45359237 kg 


1 lb(force) = 4.448222 N 


1 kip = 4.448222 kN 
1 psi = 6.894757 kPa 
1 ksi = 6.894757 MPa 


iinaki ps =— 05) 1 29648eN.m 
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SI — Imperial 


1 mm = 0.0394 in 
(imine 5 Ue ene 
(kG == 82h 20pslDmass) 
1 No=302 2258 lb force) 
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